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MODE  COUPLING  IN  THE  MODIFIED  BETATRON 


I.  Introduction 

For  small  displacements  about  a  planar  reference  orbit,  particle  motion 
transverse  to  the  toroidal  field  in  the  modified  betatron  may  be  represented  as  a 
linear  superposition  of  two  eigenmodes  of  motion:  a  "fast"  mode,  corresponding 
to  gyration  about  the  toroidal  field  lines  and  a  "slow”  mode,  corresponding  to  an 
drift  motion,  where  here  the  force  $  is  due  to  the  weak  focusing  betatron 
fields,  space  charge  forces,  and  induced  (wall  image)  fields.  The  linear  theory 
of  orbits  in  the  modified  betatron  has  been  worked  out  in  some  detail*-^  and  will 
only  be  reviewed  as  needed  here.  In  the  present  paper  we  will  mainly  discuss  the 
effect  of  quadratic  non-linearities  on  the  motion. 

Non-linear  terms  in  the  equations  of  motion  become  important  to  consider 
if:  (1)  displacements  from  the  reference  orbit  become  large,  due  say  to  the 
method  of  injection  used  or  due  to  the  operation  of  an  instability  of  some  kind, 
(2)  strong  non-linearities  (e.g.  large  values  of  3n/3r)  are  present  in  the  magnet 
design,  or  (3)  the  non-linear  term  itself  contains  a  resonant  part.  In  the 
following  we  will  illustrate  two  effects  of  quadratic  non-linearities  on  single 
particle  motion,  viz.  the  amplitude  dependence  of  the  betatron  frequencies  and 
the  exchange  of  energy  between  the  oscillation  modes  under  certain  conditions. 
These  conditions  turn  out  to  be  analogous  to  the  so-called  Walkinshaw  resonance** 
in  accelerators  without  a  toroidal  magnetic  field^-^.  We  will  limit  ourselves 
here  to  consideration  of  single  particle  motion  only,  neglecting  the  effects  of 
self  fields;  the  treatment  here  then  will  only  be  valid  for  fairly  large  values 
of  y  in  high  current  devices,  such  that  v/y  «  1*  where  v  is  Budker's  parameter. 

Four  sections  follow.  In  the  first  we  introduce  some  notation  and  sketch 
the  derivation  of  the  equations  of  motion  to  second  order  in  displacements  from 
and  transverse  velocities  about  the  reference  orbit,  taken  to  be  a  circle  in  the 

symmetry  plane.  In  the  second  section  the  equations  of  motion  are  solved 
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perturbatively  and  a  condition  for  the  generalized  Walkinshaw  resonance  is 
obtained.  Under  this  condition  we  study  the  behavior  of  the  betatron 
oscillations ,  giving  a  numerical  example  as  illustration.  As  an  interesting 
result  we  find  a  particular  value  of  field  index  gradient  for  which  the  resonance 
ceases  to  have  any  effect  on  particle  motion. 

In  these  first  two  sections  the  discussion  will  assume  that  all  applied 
fields  are  azimuthally  symmetric.  In  the  presence  of  field  perturbations  other 
orbital  resonances  may  occur  and  it  is  interesting  to  ask  whether  the  amplitude 
dependence  of  the  betatron  frequencies  induced  by  the  non-linearities  are 
sufficient  to  keep  the  oscillation  amplitudes  at  finite,  but  tolerably  small 
values.  Though  in  general  this  is  a  difficult  question  we  will  discuss  a 
special,  simple  case  in  the  third  section  of  this  paper  in  which  the  Walkinshaw 
resonance  coincides  with  both  an  integer  and  half-integer  orbital  resonance. 

A  final  section  summarizes  these  results  and  states  some  conclusions  and 


conjectures . 


II.  The  Equations  of  Motion 

The  geometry  of  the  modified  betatron  is  shown  in  Fig.  1.  We  employ 
standard  (r,e,z)  cylindrical  coordinates.  The  exact  equations  of  motion,  using 
9  in  favor  of  time  for  our  independent  variable,  are  written 


[r'(r'2  +  r2  +  z'2)“1/2]  -  r(r'2  +  r2  +  z'2)~1/2  +  x"1  (rBg  -  z'BQ)(l) 


|g-  [z'(r'2  +  r2  +  z'2f1/2]  -  X_1  (r'Bfl  -  rBr) 


(2) 


2 

where  S  is  any  suitable  function  of  position  (r,0,z),  X  5  “me  By/e,  e  and  m  are 
the  magnitudes  of  the  electron  charge  (e  >  o)  and  mass,  3  and  y  are  the  usual 
relativistic  factors,  c  is  the  speed  of  light  and  a  prime  (.')  denotes  d/d9. 

We  shall  assume  that  Bf  vanishes  on  the  plane  z  *  0  and  take  all  fields  to 
be  independent  of  9.  (The  assumption  of  azimuthal  symmetry  will  be  relaxed  in 
Section  IV,  below.)  We  take  the  equilibrium  orbit  of  a  particle  of  relativistic 
factor  yQ  at  r  •  rQ,  z  -  0,  so 


X  -  -  r  B  (r  ,o) . 
o  z  o 


(3) 


Let  us  now  define  the  normalized  coordinates  x  =  (r  -  rQ)/rQ  and  y  =  z/rQ* 
The  vector  potential  is  given  correctly  to  third  order  by 


r  B  [1 
o  zo1 


+  ^x2+yy2  -  -J  xy2  +  J  (a  +  n2  -  3) 


x3] 


(4) 


where  B  ,  n,  and  n-  are  constants.  The  corresponding  fields  are 
zo  ■ 


B  -  -  B  [n  -  n,x]y 
r  zo  2 


(5) 
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from  which  n.  is  identified  as  the  second  radial  logarithmic  derivative  of  B  . 

4m  Z 

The  toroidal  field  is  assumed  to  be  given  by 

Ba  ■  », :  V1-*  +  '  •  •>  <7) 

where  B„  is  the  value  of  the  toroidal  field  at  the  reference  orbit,  x  *  y  *  0. 

Using  the  fields  (5,6,7)  in  the  equations  of  motion  (1,2)  and  keeping  terms 
only  of  quadratic  order  gives  the  coupled  equations: 

x"  +  (l-n)x  -  by"  +  (2n-l-  ~-)xz  -  +  j  (x"  -  y" Z)  (8) 

y"  +  ny  -  -bx"  -  (2n-n2)xy  +  x"y"  (9) 

where  b  =  B.  /B  .  These  equations,  (8)  and  (9),  are  our  starting  points.  In 

U  O  ZO 

the  following  section  we  examine  the  behavior  of  an  approximate  solution  to  (8) 
and  (9)  for  various  values  of  n,  n»,  and  b. 


III.  Perturbative  Solution  of  Equations  of  Motion 

In  general  the  quadratic  terns  in  (8,9)  will  be  small  so  we  attempt  to  treat 
the  equations  perturbatively.  Neglecting  the  non-linear  terms  altogether  one  has 
the  solution  to  the  linear  equations: 


where  and  Ag  are  complex  numbers  depending  on  particle  initial  conditions  and 
the  frequencies  are  given  by 


■b  +1  -  [ (b  +1)  -  4n(l-n) 

2 


1/2  1/2 


The  subscripts  f  and  s  are  used  here  and  below  to  label  the  amplitudes  and 
frequencies  of  the  fast  and  slow  oscillation  modes.  We  will  assume  that  the 
linear  motion  is  stable,  that  is  n(l-n)  >  0. 

We  may  calculate  the  correction  to  (10)  due  to  the  non-linear  terms  by 
Inserting  (10)  in  (8,9)  and  resolving.  The  resulting  equations  will  be 
Inhomogeneous  with  various  "driving"  terms  at  the  frequencies 
2Vf ,  2vg,  0,  and  +  vg.  Consequently,  the  non-linear  correction  to  (10)  will 
remain  small  unless  it  happens  that 


the  condition  for  which,  from  (11),  being 


(12) 
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b2  -  |  [n(l-n)]1/2-l. 


(13) 


In  the  absence  of  a  toroidal  field,  (13)  is  satisfied  for  n  *  0.2  or  0.8  which  we 

identify  as  the  Walkinshaw  resonance^,  the  consequences  of  which  were  first 

12 

observed  in  cyclotrons  •  We  proceed  to  examine  particle  behavior  on  this 
resonance  in  the  modified  betatron. 

On  resonance  conventional  perturbation  theory  fails  and  one  must  resort  to 

some  other  method.  A  multiple  "time"  scale  analysis  of  the  problem  gives  a 

solution  of  the  form  (10)  in  which  the  complex  amplitudes  AF  and  A  are  no  longer 

I  s 

strictly  constant  but  vary  slowly  with  9;  they  are  found  to  obey  the  equations 


dA  * 

d r  "  iri(n,n2)AfAs  (14) 

^f  2 

d eT  *  ir2(n’n2>V  (15) 


where  2  are  two  real  valued  functions  of  the  field  indices  n  and  n2, 


2  2 
v  -  n  ?  2v  +  n 

l-TT  [2n  +  6vgZ  +  n2(  -  )] 

6v  v  -  n 

s  s 

2 

v  -  n 

-  -4  (-— - )  r2 

vf  -n 


(16) 


and  where,  on  resonance,  we  have 


v  2  -  4v  2  -  2(n(l-n))1/2. 
t  s 


(17) 


(An  asterisk  denotes  complex  conjugate  in  (14).) 

The  question  of  orbital  stab  lity  on  resonance  is  thus  reduced  to  the 


question  of  the  behavior  of  the  mode  amplitudes  in  (14)  and  (15). 

The  equations  (14)  and  (15)  may  be  completely  solved  in  a  straightforward 
manner;  the  solution  is  obtained  and  discussed  in  the  Appendix.  To  settle  the 
stability  issue,  however,  it  is  sufficient  to  note  that  there  is  a  simple 
integral  of  motion 

7  ri  |Af |2  +  7  r2  lAs|2  5  D/ri*  (18) 

and  consequently  the  motion  is  necessarily  bounded  if  r  5  Tj  >.  ®  which,  in 

fact  is  true  for  all  n  and  nj,  as  follows  from  (16).  On  this  resonance  energy  is 

simply  exchanged  back  and  forth  between  the  fast  and  slow  modes  of  motion. 

Though  we  have  argued  that  particle  motion  is  bounded  on  this  resonance  we 

have  not  in  fact  specified  a  bound  or  showed  that  the  bound  is  acceptable,  in 

terms  of  some  machine  aperture.  One  might  conjecture,  from  (18),  that  if  one  of 

T1  2  were  significantly  larger  than  the  other  then  transfer  of  energy  from  the 

more  "stiff"  (larger  r  coefficient)  mode  to  the  less  "stiff”  mode  would  result  in 

increasing  particle  oscillation  amplitude.  Hence  one  would  be  concerned  if,  from 
2  2 

(16),  either  v  -  n  or  v.  -  n.  It  follows  from  (17)  and  (18),  however,  that 
s  r 

this  can  occur  only  for  b  ■  0,  n  ■  .2  or  .8.  If  n  is  chosen  so  that  b  is 
0(1)  when  the  resonance  is  crossed  then  and  are  of  the  same  order  of 
magnitude  and  one  expects  this  resonance  to  be  quite  harmless.  For  specific 
initial  conditions  it  is  possible  to  find  a  bound  by  calculating  the  turning 
point  of  a  certain  particle-in-a-well  problem,  as  shown  in  the  Appendix. 

Curiously,  one  can  render  this  resonance  completely  inoperative  for  any 
particular  n  by  choosing  ^  so  that  and  ^  both  vanish.  From  (16)  and  (17) 
this  value  is  found  to  be 


Choosing  this  value  ensures  that  the  mode  amplitudes  remain  constant  when  passing 
through  the  "exchange”  resonance. 

We  proceed  to  illustrate  some  of  these  results  using  a  simple  single 
particle  numerical  orbit  integration.  The  algorithm  includes  the  fields  (5-7) 
but  does  not  use  an  expansion  of  the  force  or  acceleration.  Figures  2  and  3  show 
the  solutions  to  equations  (14-15)  and  (1-2)  respectively  for  the  case  n  *  0.5, 
n,,  *  n^  =  0.625,  b  *  0.5.  The  mode  amplitudes  are  strictly  constant,  no  exchange 

occurs,  and  the  particle  orbit  projection  retraces  itself  in  a  stable  manner  over 

and  over  again.  We  contrast  this  case  with  that  illustrated  in  Figures  4  and  5 
for  which  the  parameters  are  again  n  *  0.5  and  b  ■  0.5  but  now  with 
nj  11  0.  Now  the  mode  amplitudes  oscillate;  one  rises  while  the  other  falls  in 
order  to  conserve  D  (Eq  (18)).  The  oscillation  period,  from  Equation  (16)  in  the 

Appendix,  is  46.6  major  periods  for  the  particular  initial  conditions  chosen. 

The  particle  orbit  projection  now  simply  (and  harmlessly)  rotates  slowly  counter¬ 
clockwise. 

Our  conclusion  is  that  in  the  case  of  azimuthally  symmetric  fields  the 
generalized  Walkinshaw  resonance  is  quite  harmless  in  the  modified  betatron.  The 
exchange  of  energy  between  fast  and  slow  modes  is  expected  to  cause  no  major 
changes  in  the  beam  dimensions.  When  azimuthal  field  variations  are  present, 
however,  the  situation  changes  dramatically  due  to  a  coincidence  described  in  the 


next  section. 


IV.  A  Triple  Coincidence  Resonance 

When  n  -  1/2  (the  case  illustrated  in  Figures  2-5)  the  values  of  the  tunes 
vf  and  v  at  the  exchange  resonance  (13)  are  v.  ■  1,  v  »  1/2;  therefore  in  the 
presence  of  field  and  focusing  errors  the  generalized  Walkinshaw  resonance 
coincides  with  an  integer  and  a  half  integer  orbital  resonance.  This  triple 
coincidence  allows  us  to  study  in  detail  in  this  special  case  the  effect  of  mode 
coupling  (and  the  amplitude  dependence  of  the  betatron  frequencies)  on  the  orbit 
at  an  integer  and  half  integer  resonance.  Though  the  restriction  to  n  *  1/2  is 
necessary  for  there  to  be  a  true  coincidence,  if  n  is  near  but  not  exactly  1/2, 
the  three  resonances  will  be  "nearby"  and  will  occur  nearly  simultaneously  and 
the  analysis  below  should  still  hold  in  an  approximate  way. 

At  the  triple  coincidence  resonance  the  mode  evolution  equations  become 

dA 
_ s 

dQ 

^•f  2  * 

d r  -  ir2As  +  e2  +  e3Af  *  <21> 

The  e's  in  (20)  and  (21)  are  complex  constants  proportional  to  certain 
Fourier  coefficients  in  expansions  of  the  fields  and  their  gradients; 
specifically,  is  due  to  an  1  *  1  term  in  the  field  gradient,  leading  to  a  half 
integer  resonance,  is  due  to  an  l  *  1  term  in  the  field,  leading  to  an  integer 
resonance,  and  *s  due  to  an  £  *  2  term  in  the  field  gradient,  leading  to  a  "2- 
halves"  integer  resonance.  Were  they  present  alone  (i.e.  with  no  non-linearity) 
in  (20)  and  (21)  the  field  imperfection  terms  are  observed  to  lead  to  the  usual 
linear  and/or  exponential  growth  characteristic  of  integer  and/or  half  integer 
resonances.  In  the  prerjnce  of  mode  coupling  the  situation  is  much  less  clear. 
Since,  as  they  stand,  equations  (20)  and  (21)  cannot  be  solved  analytically  we 


lrlVl  +  £1A» 


(20) 
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muse  in  general  resort  to  numerical  Integration.  First  let  us  comment  on  what  we 
might  expect  to  see  in  the  solution. 

The  field  Imperfection  terms  in  (20)  and  (21)  act  roughly  speaking  as  source 
terms,  pumping  energy  from  longitudinal  to  transverse  motion.  If  this  energy 
flow  continues,  and  if  there  is  no  mechanism  to  return  this  energy  to 
longitudinal  motion  the  result  is  disastrous  —  a  linear  orbital  resonance.  Non- 
linearities,  however,  can  shift  the  betatron  frequency  of  the  resonant  particle 
off  resonance  for  some  finite  amplitude  of  betatron  oscillation.  (The  quantity 

A 

(n^  “  n2^  *s  Presumably  a  measure  of  the  frequency  shift  induced  by  a  given 
amplitude  oscillation.)  In  practice  though  one  can  not  say  a  priori  how  much 
frequency  shift  will  be  sufficient  to  terminate  the  growth  of  the  resonant  mode. 

A  numerical  study  seems  to  be  essential. 

For  a  numerical  example  we  will  examine  the  effect  of  the  non-linear  terms 
in  (20)  and  (21)  on  an  integer  resonance,  that  is,  in  the  following  we  shall  take 
■  0.  Cases  have  been  examined  numerically  for  various  other  combinations 
of  values  for  the  e's  with  no  major  differences  appearing  in  the  results. 

In  Figures  6  and  7  we  illustrate  the  mode  amplitudes  and  orbit  projection 

A 

for  a  pure  integer  resonance  with  no  mode  coupling  (n£  =  m  0.625, 

*  *005).  The  (resonant)  fast  mode  amplitude  grows  linearly  without  limit;  the 
(decoupled)  slow  mode  stays  at  a  fixed,  small  value.  The  particle  orbit  size 
(Figure  7)  consequently  grows  continuously. 

Turning  on  the  mode  coupling  changes  the  behavior  of  the  mode  amplitudes 
dramatically  but  has  little  apparent  effect  on  the  par.icle  orbit,  which  still 
appears  to  grow  to  intolerable  size.  This  is  illustrated  in  Figures  8  and  9 
where  we  see  that  the  mode  amplitudes  grow  to  a  certain  size  and  then  turn  over  - 
-  presumably  a  reflection  of  the  detuning  of  the  resonance  due  to  the  frequency 
shift.  The  "turnover",  however,  is  at  extremely  large  amplitudes  (Recall  that 


the  mode  amplitudes  are  normalized  to  the  radius  o£  the  device,  rQ.),  therefore 
the  particle  motion  appears  to  be  relatively  unaffected,  practically  speaking,  by 
changing  n^  from  5/8  to  0.  This  result  does  suggest  though  that  by 

A 

increasing  1^  ”  n2^  we  reduce  the  resonant  response  to  a  tolerable  value. 

Figures  10-11,  12-13,  and  14-15  show  our  results  for  ”  -1.,  -4.,  and 

* 

-10.  respectively.  We  see  that  as  “  n2 ^  *s  increased  the  node  turnover 
amplitude  is  reduced  and  the  particle  orbit  becomes  somewhat  more  compact, 
staying  within  the  plot  boundaries  for  significantly  longer  times.  (Even  so,  the 

a 

transverse  orbit  size  is  rather  large,  even  for  the  largest  values  of  -  n-l 
we  have  tried.) 

These  results  suggest  that  to  stabilize  the  l  m  1  resonance  a  significant 

non-linearity  (large  value  of  n^)  could  be  intentionally  introduced  in  the 

betatron  field.  One  must  be  careful  in  drawing  this  conclusion,  however,  because 

such  a  non-linearity  has  well  known  adverse  effects,  among  them  a  sensitivity  of 

the  behavior  of  the  orbit  to  initial  conditions;  that  is,  only  some  special  class 

of  particles  may  be  confined  while  others  are  lost.  Also,  if  nj,  which  is 

effectively  the  radial  derivative  of  n,  is  very  large,  it  then  becomes  difficult 

to  keep  n  itself  within  the  stable  range  0  £  n  _<  1  everywhere  within  the 

aperture.  Consequently  we  conclude  that,  as  a  practical  matter,  it  is  best  not 

to  rely  on  non-linearities  to  stabilize  the  l  *  1  resonance  in  the  modified 

betatron  and  to  design  the  machine  with  a  flat  radial  index  profile.  Avoidance 

of  this  resonance  as  well  as  other  low  order  resonances  -  -  which  then  becomes 

the  only  reasonable  experimental  alternative  -  -  is  possible  in  principle  by 

accelerating  with  constant  b  (i.e.  B.  «  B  ),  thereby  keeping  the  tunes  fixed^ 

6o  zo 

-  -  except  for  the  tune  shift  due  to  space  charge,  which  affects  the  fast  mode 
tune  only  very  slightly  for  large  B  ;  the  slow  mode  tune  can  generally  be  chosen 
to  be  very  small  (~  .2  -  .3)  for  all  time. 
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V.  Summary  and  Conclusions 

We  have  examined  Che  effects  of  mode  coupling  on  single  parCicle  orbits  in 
the  modified  betatron.  We  find  that  a  generalization  of  the  Walkinshaw 
(exchange)  resonance  can  occur  for  any  value  of  field  index  in  the  range 
.2  n  _<  .8  but  that  its  effect  on  particle  orbits  in  general  is  quite  modest 
(Figures  2-5)  and  may  be  rendered  completely  ineffective  by  a  special  choice  of 
field  index  gradient  (19). 

When  n  is  near  1/2  the  exchange  resonance  coincides  with  both  Integer  and 
half  Integer  resonances.  An  examination  of  orbit  behavior  at  this  triple 
coincidence  shows  that,  as  a  practical  matter,  the  amplitude  dependent  frequency 
shift  in  the  betatron  oscillation  due  to  mode  coupling  is  not  sufficient  to 
stabilize  the  l  ■  1  integer  resonance  (though  presumably,  as  in  the  case  of 
accelerators  not  employing  toroidal  fields,  higher  order  resonances  will  be 

O 

subject  to  non-linear  stabilization  ).  This  fact  makes  it  advisable  to  allow,  in 

the  design  of  an  experiment,  for  acceleration  with  constant  or  nearly  constant 

ratio  B  / B  thereby  holding  the  tunes  approximately  fixed  in  time. 

“O  zo 
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Appendix 

In  this  Appendix  we  discuss  the  solution  to  Equations  (14)  and  (15)  in  the 
text.  Writing 


Af  -  afexp(i<j>f)  Ag-  agexp(i<j>g) 


where  a*  -  and  a.  are  real  we  find  in  a  straightforward  way  an  equation 
1 » s  i ,  s 

.  12 
for  P  =  2  a3  : 


where 


\  P'2  +  V(p)  -  0 


3  2  12 

V(p)  -  4fp  -  4Dp  +|C 


r  "  ri  r2 


C  »  r,  a, a  cos(A_-  2<b  )  *  constant 

1  f  s  rf  Ys 


and  the  constant  D  is  defined  in  the  text,  Eq.  (18).  The  other  quantities  are 
given  in  terms  of  p  by 


a8  -  (2p ) 


“  C/(2p) 


Af  -  -g—-  (C  -  ip')  exp(2i*g). 
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These  expressions  hold  for  p  >  0.  If  p  ■  0  at  0  ■  0,  say,  then,  from  (14) 


and  (15)  A  =  0  and  A.  remains  fixed  for  all  9  >  0.  The  modal  frequency  shifts 
are  given  directly  by  (A4)  for  the  slow  mode  and  may  be  obtained  from  (A5)  for 
the  fast  mode. 

It  may  be  shown  that  V(p)  has  one  negative  and  two  positive  roots.  Denoting 
these  by  p^  ^  3  with  p^>P2>0>p^we  find  that  the  exchange  period  (period  of 
p(0))  is  given  by 


1/2 


T(Pl  -  p3)- 


K(o) 


(A6) 


where  m  *  (p^  -  p3)/(p^  -  p^)  and  where  K  is  the  usual  complete  elliptic 
integral.1^ 

The  special  cases  C  ■  0  and  r  ■  0  lead  to  motion  of  infinite  period  and 


p  ■  constant,  respectively. 


UNPERTURBED 

ORBIT 


Che  modified  betatron 
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Fig.  2:  |A  | ,  |A  I  vs.  major  periods  (A  major  period  is  a  change  of  0  by 

2it.)  n  ■  i/2;  ^  ■  5/8.  Initial  (9  ■  0)  values  of 
-3  * 

A^“  5.7735  x  10  (1  -  i)  and  Ag  ■  2A^  correspond  to  those  of  a 

particle  initially  at  x  -  y  -  .02  with  zero  transverse  velocity.  The 
same  initial  values  are  used  in  all  subsequent  figures. 
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Fig.  4: 


|A, I,  |A  I  vs.  major  periods;  n  -  1/2; 


-6 


6 


10 


z  vs. 
mode. 


*  *  ■  ■  1  1  ■  »  ■  ■ 

-2.  2. 

r-r0  (cm) 


r  -  r  ;  n  -  1/2;  n,  -  5/8;  SB  -  0.5  gauss  in  l  -  1 
o  I 

Particle  leaves  plot  area  after  ~  9  major  periods 


i 

Fourier 
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I 


•* 


I 


1  1  !  1  1  1 

i  i  1  i 

1  1  1  1  ! 

■  i  I  i 

j _ i _ 

BE 

-6. 

-2. 

2. 

6. 

r-r0  (cm) 


Fig.  13: 


2  vs.  r  -  r  ;  n  •  1/2;  n.  ■  *4;  SB  •  0.5  gauss, 
o  i 

plot  area  for  50  major  periods. 


Particle  remains  in 


% 
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